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A system of trapped ions under the action of off-resonant standing-waves can be used to simulate 
a variety of quantum spin models. In this work, we describe theoretically quantum phases that can 
be observed in the simplest realization of this idea: quantum Ising and XY models. Our numerical 
calculations with the Density Matrix Renormalization Group method show that experiments with 
ion traps should allow one to access general properties of quantum critical systems. On the other 
hand, ion trap quantum spin models show a few novel features due to the peculiarities of induced 
effective spin-spin interactions which lead to interesting effects like long-range quantum correlations 
and the coexistence of different spin phases. 



I. INTRODUCTION 

The goal of quantum many-body physics is the under- 
standing of properties of materials ranging from quantum 
magnets to high-Tc superconductors. For this task, the- 
orists have developed a variety of simplified models which 
allow us to describe the rich phenomenology that is ob- 
served in those systems. In the last years, progress in 
the manipulation of matter at the microscopic scale has 
changed dramatically this situation. Atomic experimen- 
tal systems such as cold atoms in optical lattices offer us 
the possibility to tune interactions and engineer at will 
quantum phases. In this way quantum interacting mod- 
els can be studied in a controlled way and some of the 
limitations of solid-state set-ups are overcome. A suc- 
cessful step in this direction was given recently with the 
observation [| of the superfluid-Mott insulator transi- 
tion with bosonic atoms in an optical lattice. Recently, 
cold atoms in optical lattices have also been proposed for 
the simulation of certain quantum spin models 

On the other hand trapped ions are a clean experimen- 
tal system where quantum optics offers us accurate tech- 
niques for the manipulation and measurement of quan- 
tum states 3J. In the last years this experimental field 
has been mainly motivated by applications to quantum 
information processing 4], in which internal electronic 
states are used as qubits, and vibrational modes permit 
us to perform quantum gates between them 0,0- In 
fact, by using similar ideas we have shown recently that 
trapped ions can be used for the study of a rich variety 
of quantum interacting models Q, Q ■ The experimental 
requirements for the study of quantum interacting sys- 
tems with trapped ions are indeed much less stringent 
than those for quantum information tasks. 

In this paper we pursue the ideas proposed in Ref. 
for the realization of a quantum simulator 0,0] of quan- 
tum magnetism with trapped ions. We focus on the pos- 
sibility of studying antiferromagnetic Ising and XY spin 
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chains in linear ion traps, where quantum phase transi- 
tions can be induced and explored We show that 
experiments with trapped ions can access the physics of 
magnetic quantum phases with an accuracy that is not 
possible in other experimental set-ups, and allow us to 
test general properties of quantum systems at, or near, 
criticality. On the other hand, the effective spin Hamil- 
tonians that can be engineered with trapped ions show 
a few remarkable features, the most important one be- 
ing the fact that spin-spin interactions follow power-law 
decays, something that induces long-range quantum cor- 
relations between distant ions, which are absent in the 
usual nearest-neighbor models. 

Our paper is organized as follows. In the first section 
wc describe in detail our proposal, taking as a starting 
point the system formed by a set of trapped ions cou- 
pled to an off-resonant standing wave, and show that it 
realizes a system of effective interacting spins. We fo- 
cus on the cases of Ising and XY interactions. In the 
second section we study the many-body problem posed 
by the quantum Ising model in ion traps. For this pur- 
pose, we make use of a quasi-exact numerical method, 
the so called Density Matrix Renormalization Group ^2] j 
and introduce also an analytical, approximate solution 
in terms of Holstein-Primakoff bosons. The third sec- 
tion deals with the study of the XY model in ion traps. 
Finally, we present our conclusions. 



II. EFFECTIVE SPIN QUANTUM SYSTEMS 

In this section we review the theoretical description of a 
set of trapped ions under the action of off-resonant stand- 
ing waves presented in , and show that their quantum 
dynamics follows approximately that of quantum inter- 
acting spins. We also discuss carefully the effects of resid- 
ual spin-phonon couplings. 

Our proposal relies on an always-on coupling between 
internal states and vibrational modes, in such a way that 
phonons transmit an interaction between different ions. 
This paper deals, thus, with the theoretical description of 
a system of effective spins (internal states) coupled to a 
set of vibrational modes. The corresponding Hamiltonian 
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has three contributions: 



H = H m + H v + H f . 



(1) 



H m describes the local dynamics of internal states, which 
are a set of independent effective spins under the action 
of effective global magnetic fields in each direction: 
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j=l,...,N 
a.—x.y,z 



B a af. 



(2) 



Note that B z is the energy of the internal state, and B x , 
B v can be implemented by the action of lasers resonant 
with the internal transition. In Eq. IjTJl. H v is the vibra- 
tional Hamiltonian, and Hf is the internal state-phonon 
coupling. Both terms are described in the following two 
subsections. 



A. Vibrational modes of ion chains 

We consider ID systems of trapped ions, whose phys- 
ical implementation corresponds to Coulomb chains in 
linear Paul traps or linear arrays of ion microtraps. Let 
us assign z to the axis of the ion chain, and x, y to 
the radial directions. The potential experienced by the 
ions is determined by the trapping frequencies, Lu a (in 
the following the index a always runs over the spatial 
directions, a — x,y, z), and the Coulomb repulsion: 
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In the harmonic approximation we express V as a func- 
tion of the displacements around the equilibrium posi- 
tions (q" = x" — up to quadratic terms: 
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The corresponding Hamiltonian can be diagonalized in 
terms of collective modes (phonons): 

1 p a 
H v = 2 m E K-tilitfi +E 2^ = E ^^.A"! 5 ) 



where pf are the momenta corresponding to qf. Local 
coordinates can be expressed in terms of creation and 
annihilation of collective vibrational modes: 
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where c x . y = 1, c z = —2. Matrices A4 a in llOj l diagonal- 
ize the vibrational hamiltonian: E< , Mf n lC"jAi? m — 



J a.n u rL,m- 



The characteristics of the vibrational modes are gov- 
erned by the following parameters, which quantify the 
relative value of Coulomb interaction and trapping po- 
tentials: 
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where do is the mean distance between ions. If (3 a <C 1, 
then phonons are close to be localized at each ion (stiff 
limit); on the contrary j3 a S> I (soft limit) implies that 
Coulomb repulsion prevails over the trapping potential 
and, thus, phonons have a strong collective character that 
results in the ability to mediate interactions with a range 
of the order of the ion trap ■ 

We summarize below the description of vibrational 
modes in different trapping conditions: 

(i) Coulomb chains (Paul traps). In this experimen- 
tal set-up the equilibrium positions in the axial (z) 
direction are determined by the Coulomb interac- 
tion and the axial trapping, and [3 Z depends only 
on the number of ions N . This dependence can be 
estimated in the limit of many ions (see |13| ) : 
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(3 Z ^> 1 and thus axial modes in Coulomb crys- 
tals are always in the soft limit. On the other 
hand, /3 x ,s/> can be reduced at will because one can 
increase the axial trapping frequencies u) XtV while 
leaving constant the mean distance between ions, 
d . Indeed, condition f3 xy < 1 has to be fulfilled 
to make the Coulomb chain stable against zig-zag 
instabilities |l3[ . 

(ii) Linear arrays of ion traps. It is worth consider- 
ing this case here due to the experimental effort 
currently being devoted to the fabrication of linear 
arrays of ion traps where equilibrium positions of 
the ions are chosen by individual confinement po- 
tentials for each ion |14|. In this case, all trapping 
frequencies u> a can be chosen at will, and in partic- 
ular condition /3 a < 1 can be also fulfilled opposite 
to the case of standard Paul traps. 



B. Internal state conditional forces 

Internal states are coupled to the motion by placing 
the ions in an off-resonant standing-wave, such that they 
experience a state-dependent a. c. -Stark shift. Several 
configurations are possible, in which internal states are 
coupled to vibrational modes that are transverse or longi- 
tudinal with respect to the trap axis. The characteristics 
of the resulting effective spin-spin interaction depends 
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on the choice of the directions of the laser beams: in 
particular, axial vibrational modes mediate long-range 
interactions, with a range of the order of the ion chain, 
and radial modes mediate shorter range spin-spin inter- 
actions, with a power-law decay Ji.j oc — j| 3 . Due 
to the fact that short-range spin Hamiltonians contain 
a rich quantum critical phenomenology, we focus in this 
work on this last situation. 

Let us consider the following coupling between effective 
spins and radial motion: 
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(9) 



\l) a is the eigenstate of a a with eigenvalue 1. The rea- 
son for this choice of couplings is that conditional forces 
are most easily implemented with the z component of the 
internal states (one only needs a single standing-wave). 
Thus, it is advantageous to couple a z to one of the radial 
directions. The first term in Hf is the one which corre- 
sponds to the pushing gate presented in |l5j |. while the 
second one can be implemented by using two additional 
standing-waves Q- 

The coupling of the internal states to the motion can be 
written as an effective spin-phonon coupling by express- 
ing the ions' coordinates in terms of collective modes: 
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In Eq. UlUfl w e have introduced the notation a x = cr z , 
o-y = a y , that is, a a is the component coupled to the 
motion in the direction a. 



C. Effective spin— spin interactions 

Under certain conditions, a set of spins coupled to a 
common bath of vibrational modes is well described by 
an effective spin interacting Hamiltonian, something that 
is apparent if one makes use of the following canonical 
transformation [7lll6|: 
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where rjf n are the displacements of the modes in units of 
the ground state size. 

In the new basis the Hamiltonian JJJ includes an effec- 
tive spin-spin interaction: 

e- s He s = H v + \ £ J^af a«+]T B' a a?+H E , (12) 



where 



J h3 ~ m, ,2 Jvl i,n Jvl 3,n — z 'li,n'lj,n nw a,n- 
n ' n 

(13) 

He is a residual spin-phonon coupling, whose explicit 
form is given below. From (|13fl and the definition of A4 a 
one can easily deduce that in the limit [3 a -C 1: 



= 2(3 a r) 2 a t>u a , (14) 



where z 1 i are the ions' positions in units of the average 
distance, do. We have introduced 77, which characterizes 
the displacement of the phonons due to the presence of 
the state-dependent force: 



a = F a y / h/2muj a /huj a . 



(15) 



Note that in 114(1 the effective magnetic fields receive, 
after the canonical transformation, a contribution from 
the pushing forces B' a = B a + l^/(muj). The extra 
term in B' a does not depend on the ion's site, and thus 
can be considered as an overall correction to the global 
effective magnetic fields. 

Radial modes in a chain of trapped ions allow us to 
implement two types of spin models. If Hf acts only on 
one of the directions of motion, say x, we get an Ising 
spin-spin interaction: 
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The residual spin-phonon couplings are given, to lowest 
order in 77 by: 
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On the other hand, if we apply two conditional forces 
in both radial directions, we get a Hamiltonian which 
couples two components of the effective spins. 

i,j i 

A magnetic field B x is included such that, in a rotated 
basis a z — > a x , a x — > — g z , and l|18|l corresponds to the 
XY model in its usual basis. The residual spin-phonon 
couplings are given, to lowest order in r\ a by: 
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( a a,n + a a,n)( a ' < a ' m "a' ,ro) 



.(19) 



Hg Y accounts for the interference between the two differ- 
ent conditional forces, which can be avoided by choosing 
different values for the two radial trapping frequencies, 
as we explain in the following subsection. 
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D. Decoherence induced by vibrational modes 

The residual spin-phonon couplings (|T71 IT§|) are a 

source of decoherence which entangle internal states with 
phonons. Besides that, internal states are prepared and 
measured in a different basis than the one correspond- 
ing to effective interacting spins. Both effects have to be 
evaluated to study how the quantum dynamics of inter- 
nal states departs from the ideal quantum spin Hamilto- 
nian. The problem of solving the quantum dynamics of a 
system of interacting spins coupled to phonons is indeed 
quite complicated. We are interested here in estimating 
the error induced by these couplings, so that we will make 
use of perturbation theory and several approximations. 

Let us consider that the system is initially in a prod- 
uct state of a thermal phonon density matrix, p p h, and 
a given internal pure state, pi = \ipi)(ipi\. \ipi) evolves 
to \ipf) under the spin Hamiltonian iJg, thus, \ipf) rep- 
resents the ideal simulated spin state. We define the fi- 
delity in the quantum simulation by the overlap between 
the state of the system after the evolution with the whole 
Hamiltonian, and \ipf): 

Ftyi) = (^ f \TT ph {e-' lHt ^ Pl ® p ph e lHt l h }\^ } ). (20) 

Having in mind perturbative calculations, let us express 
the fidelity in the following way: 

F{<M = (21) 
(^|Tr P 4e 5 ( t )[/(t)e- 5 ( )p l ®p p ,e 5 (°)C/(i) t e- s( * ) }|V J ), 

where U(t) = e iH t/h e -i{H +H E )t/h^ is the evolution op- 
erator in the interaction representation with respect to 
He- Hq = H v + Hs is the Hamiltonian without resid- 
ual spin-phonon coupling. In Eq. 1|2I[) . as well as in the 
right-hand side of the equations below, all the operators 
evolve with Hq. 

We are particularly interested in the very practical 
question of how do measurements of internal state ob- 
servables departure from the averages in the simulated 
quantum spin models. Let us consider a few-site spin 
operator, O, and define the error in its average, Eo — 
(0(t)) - {O(t)) , with: 

<0(i)> = 

Tr{0(t)e s Wu(t)e- s ^ Pl (g>p ph e s ^U(t^e- s W}, 
(O(f)> = (*/|0|*/>. (22) 

Equations (|2f(l and (|22() provide us with a well suited 
starting point for calculating the effects of the residual 
spin-phonon coupling in a series in rj a . We show below 
that these terms can result in negligible contributions 
with the right choice of parameters. Due to the fact that 
lowest order spin-phonon couplings are different in each 
quantum spin model, we consider the two cases sepa- 
rately. 



1. Ising Model 

The spin-phonon couplings in I jl 7|) are proportional to 
the transverse magnetic field, B x . In the following we 
ignore the index a in the vibrational modes, and assume 
that it corresponds to one of the radial directions. In 
the most interesting regimes B x sa J, thus we estimate 
B x « rfuj. All terms in H E sm9 are non-resonant, such 
that, if we ignore the contributions from S in (1211 122fl . 
the only allowed transitions are virtual with probability 

The most important contribution to the error is thus 
the one that comes from the canonical transformation 
only, which can be estimated by setting U(t) = 1, ex- 
panding e s '*' in l|21|). and keeping terms up to order i] 2 . 
In the stiff limit, that is (3 <C 1, we make the additional 
approximation of considering vibrational modes in 12I|) 
as localized phonons, and get the following expression for 
the error to lowest order in ty. 

E = 1-^ = E^+E^-E^-E^°\ 
E [tiM = v 2 J2( fieiU ^ t2) + ( fi + 1 ) e ~ iu(tl ' t2) ) x 

3 

«<7|(t 1 )<7?(t 2 ))-<a|(t 1 ))<a|(t 2 ))),(23) 

where n is the mean radial phonon number. The error 
is of order Nrj 2 , and is proportional to the fluctuations 
in cr?. Note that equal-time averages give the dominant 
contribution in (|23[) because the two-time spin averages 
m 

E [o,t] 

can be neglected at long enough times. 
The scaling of E with N results from the fact that 
T in (|2I|I represents the overlap between the internal 
state of the ion chain and the effective spin state. How- 
ever, in practice, local observables like single spin aver- 
ages or two-site correlation functions are measured. We 
can show that in this case the error does not increase 
with N. Let us consider again the effect of the canonical 
transformation in Eq. (|22() . We get the following error 
in the measurement of the observable: 

p _ p [o,o] p [t,t] 9P [o,t] 
Eo — &o + ~ Ztj O ' 

E [ o lM = l{[[0(t),S(h)],S(t 2 )}). (24) 

An explicit expression can be derived from l|24|l which is 
not very enlightening. However we note again that equal 
time correlations in l|24|) are the dominant contribution, 
such that: 

E o « ^ a E( 2R + 1 )<[[°(*)'°i( i )W(*)]). (25) 

3 

where we have also approximated vibrational modes by 
localized phonons. It is clear from Eq. (|25|) that if O is an 
M-site observable, there are only a maximum of M non- 
vanishing commutators and thus E « Mrfuj. The most 
meaningful physical quantities in the study of quantum 
criticality are indeed one-site (mean values) or two-site 
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(correlation functions) averages, and i12">l) implies that 
these ones can be studied with an error that is indepen- 
dent on the number of ions. 

2. XY Model 

In this case special care has to be paid to the effect of 
the residual spin-phonon couplings. H^ Y in (|19J) is of 

order ij^u> a , that is, of the same order as itself. On 
the other hand, if lu x = u> v , then He contains resonant 
terms that couple vibrational modes in different trans- 
verse directions. Under these conditions the effect of He 
is comparable to that of the spin-spin interaction and 
the quantum simulation is ruined. 

A way out of this problem is to tune lo x ^ u> y , 
with uj x — Lu y of the order of lo x , lo v . In this case, 
there are no resonant terms in i|19|) . and the interfer- 
ence between standing-waves in different radial direc- 
tions is suppressed. The error is then of the order of 
(rj 2 uj x ^y) 2 1 '(u> x — toy) « (r] x ,r]y). Indeed, Paul traps are 
usually designed such that the radial frequencies are dif- 
ferent, with parameters that fulfill the conditions for the 
rotating terms in H§ Y to be neglected [2J- Under these 
conditions, the lowest order contribution to the error 
is, again, the one that comes from the change of basis 
(E oc r)luj a ). 

III. EFFECTIVE ISING MODEL 

We have shown that antiferromagnetic long-range 
Ising models can be realized in experiments with ion 
traps. Hamiltonian 1)16(1 is exactly solvable in the case 
of interaction between nearest-neig hbors (NN) mm 
In this case, the sign of the interaction is not relevant at 
all, because the transformation: 

U = II ^) • °*i "> Ua]U- x = (26) 

j odd 

maps the ferromagnetic into the antiferromagnetic 
model. 

The exact solution of the NN-Ising model shows that 
there exists a quantum phase transition at B x = J 01 
between a paramagnetic state (B x > J), and an anti- 
ferromagnetic phase (B x < J) characterized by the Neel 
order parameter: 

°iv^£(-l)X- ( 27 ) 

3 

We expect that the properties of the effective spin model 
in ion traps (|16|l are similar to those of the NN-Ising 
model, due to the fast decay of interactions Jjj oc 
— j\ 3 ; in particular, we expect a quantum phase tran- 
sition at a given critical value of the longitudinal mag- 
netic field, B x . In the following we present numerical cal- 
culations to describe quantitatively the quantum phases 



of Hamiltonian i|16|) . and show that critical properties 
are very similar to those of the NN-Ising model. Indeed, 
renormalization group arguments can be used to show 
that the 1/r 3 Ising model belongs to the short-range 
Ising universality class • O n the other hand, the long 
range character of the interactions turns out to induce 
intriguing effects in the spin quantum correlations which 
are explained below by means of a spin-wave model. 

This numerical problem is handled with the Density 
Matrix Renormalization Group (DMRG) method [l2fl . 
which is a quasi-exact numerical method for the study of 
ground states of interacting quantum systems in ID. The 
fact that we have further than nearest-neighbor terms in- 
creases the complexity of the algorithm by a factor of the 
number of sites, N. We keep m — 128 eigenstates of the 
reduced density matrix at each step in the DMRG algo- 
rithm and test the accuracy of our calculation by com- 
paring its results with the exact solution in the NN-Ising 
model at the critical point, where DMRG works worst. 
We have found that the relative error of the ground state 
energy in this case is limited by the machine accuracy 
SE ~ 10~ 15 . One can expect the same accuracy in cal- 
culations with the ion trap Ising models presented below 
because correlations are similar in both cases. 

Our numerical calculations describe two types of quan- 
tum Ising models: 

(i) 1/r 3 -Ising interaction. 

This case corresponds to a linear array of micro- 
traps, where equilibrium positions of the ions are 
approximately constant and fixed by the position 
of the microtrap. 

(ii) Linear trap-Ising model. 

If the trap is in the stiff limit, then effective spin- 
spin interactions decay like l/\z® — z j| 3 ; but dis- 
tances between ions in the Coulomb crystal depend 
on the position. Thus, we get an effective inho- 
mogenous quantum Ising model with interaction 
strength Ji j which depends on the position (see 
Fig. The ground state shows the coexistence 
of different phases in different locations of the ion 
trap. 

The results presented in the following subsections 
where obtained with chains of TV = 100 ions. Ji.j in the 
linear trap-Ising model was calculated with (3 X = 10~ 2 , 
however, our results do not change much with [3 X , as long 
as one considers values within the stiff limit (f3 x < 0.1). 
If we assume typical values oj x — 10 MHz, and rj 2 = 10~ 2 , 
this would correspond to interaction strength J^.i+i < 20 
kHz. 



A. Effective magnetization 

In the following we study the effective magnetization 
and its fluctuations. We will be mainly interested in (i) 
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whether global measurements are enough to character- 
ize quantum phases, and (ii) which is the effect of in- 
homogeneity and finite size in linear trap models. All 
our results are presented in energy units Jo, where Jo 
is the interaction strength between nearest-neighbors in 
the l/r 3 -Ising model, or the averaged nearest-neighbor 
interaction l/N^i Ji.i+i, m the case of the linear trap- 
Ising model. 



1. Transverse ((o~ x )) magnetization. 

The mean magnetization, m x = (1/N)(a^) (where 
<7y = ^jCr"), can be obtained from global measure- 
ments. In Fig. n ( a ) we show the evolution of m x 
with the magnetic field. The magnetization curve of 
the l/r 3 -Ising model is similar to the NN case. The 
quantum phase transition results in a discontinuity in 
d 2 m x I d(B x ) 2 (see Fig. (a)), something that allows 
us to locate the critical point B x (l/r 3 ) w 0.83, which 
lies below the critical point in the NN- Ising model 
(B X (NN) = 1). This effect can be explained in terms 
of frustration induced by terms J^j, with (i — j) even, 
which reduces the stability of antiferromagnetic order. 

On the other hand, m x in the linear trap-Ising model 
departs from the homogeneous I/r 3 -Ising case, due to 
the spatial variations in Jjj (see Fig. d(b)). The system 
shows the coexistence of different phases, depending on 
the local value of Jij, as shown in the evolution of the 
local magnetization (a x ) with B x in Fig. [3 (b). The 
local phase diagram at each ion j shows a critical point 
B x (j) w 0.83Jjj+i, which is determined by the local 
value of the interaction, but satisfies the same relation 
with Ji i+i found for the homogenous l/r 3 -Ising model 

(Fig. m'(b)). 

The local phases in a linear ion trap Ising model can 
be studied with only limited local addressing of the ions. 
For example, by measuring the average magnetization, 
m x , corresponding to the 20 central ions, we can observe 
the signature of the quantum phase transition, as shown 
in Fig. |21 (a). Thus, individual ion addressing is not 
necessary for detecting the critical point in the quantum 
simulation. 



2. Longitudinal ((o~ z )) magnetization. 

The antiferromagnetic order parameter is given by the 
staggered magnetization <|27|) . In the thermodynamic 
limit, spontaneous symmetry breaking results in a non- 
zero value of (On)- On the contrary, in finite systems, 
symmetry remains unbroken. In this case, it is convenient 
to study the squared antiferromagnetic order parameter 
(O 2 ^), which takes a value rs 1 in the antiferromagnetic 
phase. In Fig. |3 we present the evolution of (O n ) in the 
1 /r 3 -Ising case, as well as in the central region (20 ions) 
of a linear trap. 
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FIG. 1: (a) Evolution of the averaged effective longitudinal 
magnetization in 1/r 3 and linear trap Ising models (JV = 100 
ions). We also plot (crj) in the case of the central ion (j = 50) 
in a linear trap, (b) Local strength of the Ising interaction in 
a linear ion trap, and local value of the critical field. 
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FIG. 2: (a) Second derivative of m x as a function of B x in 
ion trap Ising models (N = 100 ions). Dashed line: 1/r 3 - 
Ising model. Solid line: average over the 20 central ions in a 
linear trap-Ising model. Dotted line: central (j = 50) ion in 
a linear trap-Ising model (solid and dotted line are almost on 
top of each other), (b) Contour plot of (<rj) on the plane of 
the coordinate j and the magnetic field B x . 



For measuring On it is necessary to address each ion 
individually. On the other hand, the fluctuation of the 
average longitudinal magnetization, (1/N 2 )((a^) 2 ) is an 
interesting alternative which does not require individual 
ion addressing. Antiferromagnetic order can be detected 
by means of this global observable, because longitudinal 
spin fluctuations are suppressed in the Neel ordered state, 
as shown in the case of the 1 / r 3 -Ising model (Fig. [2J • 



B. Correlation Functions 



Correlation functions Cff = (afaj) - (af)(af) char- 
acterize quantum correlations in the ground state of the 
effective spin system. Even when they play an important 
role in the description of strongly correlated systems, it 
is not possible to measure them directly in solid-state 
experimental setups [l9j . Trapped ions, on the contrary, 
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FIG. 3: (a) Fluctuation of the longitudinal effective magneti- 
zation, and (b) Neel order parameter. In both plots we con- 
sider the l/r 3 -Ising model, and the central region (20 ions) 
of the linear trap-Ising model. 
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FIG. 4: (a) Absolute value of correlations, \C 



jo do +j I ' 



50 



between 



the central ion (jo = 50) and the rest of a chain with N = 100 
ions, in the case of the l/r 3 -Ising model, (i) C®* j +j, B x = 
0.72 < B x ; (ii) C% do+j , B x = 1.32 > B x . (b) Zoom of the 
long-range tail of C zz , which follows an algebraic decay with 
an exponent a = — 3. 



offer us the possibility to measure directly equal-time 
correlation functions by means of a set of measurements 
on single ions. We show here that realizations of Ising 
models with trapped ions allow us to test directly prop- 
erties of quantum critical systems such as the algebraic 
decay of correlations at a quantum phase transition, as 
well as remarkable new effects induced by long-range in- 
teraction terms. 

We consider correlations of observables that are trans- 
verse to the order parameter, that is, Cf J correlations in 
the antiferromagnetic phase and Cfj correlations in the 
paramagnetic one, because they are the most meaningful 
in terms of the spin-wave picture to be introduced later. 
In Figs. 01 and it is shown that quantum correlations 
present two regimes, depending on the distance between 
ions, \i - j\: 



2. Long-range correlation mediated by the interaction. 

A remarkable feature in Figs. 0] and 03 is that correla- 
tion functions decay like a power-law, CfJ w l/\i — j\ 3 , 
at very long distances, so that properties of quantum cor- 
relations depart from the nearest-neighbor case (see Fig. 

H(b)). 

This effect can be qualitatively understood by consid- 
ering that entanglement between distant ions can be cre- 
ated in two ways: (i) by nearest-neighbor terms in the 
Hamiltonian, in such a way that correlations present the 
same characteristics as those of the NN-Ising model (ex- 
ponential decay), or (ii) directly by long-range terms in 
(|16|l in such a way that they decay following the power 
law of the spin-spin interaction. This hand-waving ar- 
gument will be justified in the following section by means 
of a spin-wave formalism. 



1. Universality regime. 

At intermediate distances, correlation functions both 
in 1/r 3 , and in linear trap-Ising models, show critical 
properties which are similar to those of the NN-Ising 
model: Cf" cx e~^~^^ a " when B x ^ B%, whereas 
Cff cx \i - j\- p at B%, with p = 2. On the other 
hand, correlation lengths £ aa (see Fig. 01 diverge near 
the critical point and show the dependence, (£, aa )~ 1 cx 
\B X — B x \ u , with v ~ 1. The critical properties of 
the l/r 3 -Ising model are thus the same as those of the 
nearest-neighbor model [l7j . 

We note that by measuring quantum correlations in 
the central region of the chain one can measure critical 
exponents even in the case of the linear-trap Ising model. 
Finite size effects are however more important in this case 
(see Fig. 03(b)). 



C. Spin— wave picture 

Far away from the critical point, excitations over the 
ground state can be described as spin-waves satisfying a 
harmonic Hamiltonian. This picture allows us to under- 
stand the numerical results presented above; in particular 
correlations in the presence of power-law interactions. It 
also brings ion trap spin models in connection with recent 
results on quantum correlations and criticality in bosonic 
gaussian ground states poj . 

Spin-waves are defined by the Holstein-Primakoff 
(HP) transformation |l9j| . whose particular form depends 
on the ground state. Thus we have two consider sepa- 
rately the two limits corresponding to antiferromagnetic 
or paramagnetic phases. 
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FIG. 5: (a) Absolute value of correlations, |C^™j +j|> between 
the central ion (jo = 50) and the rest of a chain with N = 100 
ions, in the case of the linear trap-Ising model, (i) C xx , B x = 
0.89 < B x - (ii) C zz , B x = 1.72 > B x . (b) Plot of C XX >J0+J for 
both 1/r 3 and linear ion trap models, exactly at the critical 
point B x = B x . 



Let us write this Hamiltonian in terms of canonical op- 
erators: 

Qi = -^$ + 6 0, 

Pi = -^(bl-bl), 

H/\2B X \ = lj2 K i,iQiQi + lj2 P ?' ( 3 °) 

3,1 i 

with Kj i = Jj,i/\B X \ +Sjj. In the following we consider 
the limit N — > oo, so that we can get analytic results. 
In this limit, Hamiltonian (|3(JI) is diagonalized by plane- 
waves, Q q = l/s/NY^jS^Qj and correlation functions 
are given by: 



Cft = 2(Q Q Q J ) = -J2< 
q 

O ■ 



"i j 



(QqQ- 



(31) 
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f2„ is the spin-wave dispersion relation: 



n 



9 = ^E^°< 



(32) 



It is illuminating to write the correlation function in the 
following way: 

os3 - W) = ^ E e ~ i9 ^ = E = 



A: 



E 



•4* 

B 1 



A l i 



(33) 



FIG. 6: Correlation lengths (^ QQ )~ 1 as a function of B x in 
both 1/r 3 and linear ion trap models, (i), (i') curves corre- 
spond to {S, xx )~ 1 , whereas (ii), (ii') correspond to (£ 22 ) -1 ). 
£ zz is not shown in the antiferromagnetic phase because CfJ 
tends to a constant value in the absence of symmetry break- 
ing. 



1. B*>J 

In this phase, HP bosons describe spin-waves excited 
over the paramagnetic ground state: 



with: 



tf-io*)/2 = 6t(i- 6 t 6j) i/2 
(of +irf)/2 = (l-6t & .)V2 6 . 



3' 



= 2b)b 3 - 1. 



(28) 



The harmonic approximation is valid if b^bj <C 1. In this 
limit the effective spin Hamiltonian takes the following 
form: 

H = \ E JjAb}+bM+bi) + B x J2( 2b % - !)■ ( 29 ) 

3,1 3 



A z — 
A i - N 



-y 



1 

27 



dq 



n 2 



(34) 



The function A? — > as j — > oo, such that for long 
distances, Eq. I|33|l implies that correlations decay fol- 
lowing the spin-spin interaction power law. Indeed, the 
decay of correlations with the same power-law than the 
interaction term, has been recently shown to be a gen- 
eral property of ground states of harmonic lattices with 
long-range interactions [2fj| . 



2. B x < J 

In this limit the ground state is close to the antiferro- 
magnetic (Neel state), so that we use the following HP 
transformation: 



2b]b 



j»3 



1, 



(of +*(-iy jf)/2 = &J(1-&&)V3 

(o?-i(-l)M)/ 2 = (1-&}M 1/2 ^ 



(35) 
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Signs, (— iy , have to be added to choose the Neel state as 
the ground state of the HP oscillators. The Hamiltonian 
is: 

H = ^E(- 1 ) J '~' J i.' { 2b h - x ) ( 26 N l) + 

3,1 

sx E( 6 l+ 6 i)' ( 36 ) 

If we neglect interactions between bosons, then we get a 
set of non-coupled harmonic oscillators. However, non- 
quadratic terms in H36|l induce correlations due to the 
presence of the transverse field B x . To see this, we solve 
first the quadratic part of the bosonic Hamiltonian: 

H =j2 2jb h+ BX H( b ]+ b i)> ( 3? ) 

3 3 

where J = — I) 3 ' - Jj,h that is, the mean longitu- 

dinal magnetic field. Hamiltonian (|37() is solved by dis- 
placing the HP bosons, which corresponds to consider 
the mean-field ground state as a reference state for the 
HP transformation: 

bj^bj-B'/pJ). (38) 

Up to quadratic terms in the displaced HP bosons, we 
get the following Hamiltonian: 

H = 2J^b]b j + 

3 

(ff) a B- 1)W " I>JA ' ( b l +b i) (^ + 6l).(3») 

3,1 

By following the same steps as in the previous case, we 
show that: 

°03 = A * + EMr'^ir 1 ^' (40) 

with A? given by the same function (|34[l of the spin-wave 
energies that diagonalize 1)3 7|) . Thus, to lowest order in 
B x /J, x-x correlations behave in a similar way as z- 
z correlations in the antiferromagnetic phase, with the 
only difference being the alternation in the sign. 

In Fig. U| we show the comparison between results ob- 
tained by means of the spin-wave picture and DMRG 
calculations in both antiferromagnetic and paramagnetic 
phases. We have checked that the agreement is good far 
away from the critical point, and gets worse when one 
approaches it because the assumptions behind the HP 
approximation are no longer valid. 

IV. EFFECTIVE XY MODEL 

We consider the isotropic case of Hamiltonian Q18JI. 
j\ x ] — j\ y ] > 0. Let us work in the rotated basis defined 




j (ion site) 

FIG. 7: (a) Absolute value of correlations | C_^,™j- + j | {jo ~ 50, 
N — 100 ions) in the case of the l/r 3 -Ising model. Points: 
Numerical results from the DMRG calculations. Solid line: 
calculation with the HP transformation. 

below Eq. ggj, with B z -> -B z : 

H s Y = \Y. + + E ( 41 ) 

i,3 i 

but will keep in mind that in experiments effective spin 
observables must be measured in the original basis. Note 
also that [H, ^ . cr|] = and hence crj is a conserved 
quantity. 

The nearest-neighbor (NN) XY model can be exactly 
solved by a Jordan- Wigner mapping to free fermions [2l| . 
Note that here we consider the antiferromagnetic model, 
which can be mapped onto the ferromagnetic one by 12011 . 
The whole region \B Z \/J < 1 is critical and CfJ follow a 
power law with critical exponent a = 1/2. (a x ) = (cr y ) 
— 0, whereas (a z ) grows as a function of B z up to the 
non-analytical point \B Z \/J = 1. One can expect that 
properties of XY models in trapped ions are similar to 
those of the NN case. 

In this section we study the two cases of 1 /r 3 -XY inter- 
actions and linear trap-Ising interactions (see the begin- 
ning of section l|IIIfl for a motivation of this distinction) 
by means of DMRG. m = 128 eigenstates of the reduced 
density matrix are kept, and comparison with the exact 
solution allows us to estimate a relative error in the cal- 
culation of the energy, SE ~ 10 -13 . We consider ion 
chains with N = 50 ions, and the same parameters for 
the ion linear trap, and energy units, considered in the 
previous section. 



A. Effective magnetization 

The most interesting single spin observable is the mag- 
netization in the z direction. Here we will follow the 
same lines and definitions explained in subsection IIII Al 
In Fig. [S] we plot the evolution of m z as a function of 
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FIG. 9: (a) Local critical field calculated from the evolution 
of (crj) in the linear trap-XY model, (b) Spatial dependence 
of {cr|} in a linear-trap XY model with different values of B z . 



FIG. 8: Phase diagram of 1/r and linear trap-XY models. 
We plot the averaged global magnetization, as well as the 
evolution of (<t|) at the center of a linear trap (N = 50 ions). 



B z . The steps in the curve are due to the finite size of 
the ion chains. The magnetization curve of the 1/r 3 - 
XY model follows the same relation as in the nearest- 
neighbor model, m z oc (B z ) 2 . We find again the same 
effect that in the Ising model, that is, the critical point 
is shifted B z (l/r 3 )(^ 0.9) < B Z {NN){= 1) due to frus- 
tration induced by long-range interaction terms. 

On the other hand, in the linear trap-XY model, m z 
departs from the homogenous 1/r 3 case, due again to the 
variations of the interaction strength along the ion trap. 
As we did in the case of the Ising model, we get a local 
phase diagram by plotting the evolution of the single site 
magnetization as a function of B z . An interpretation of 
our results in terms of local quantum phases governed 
by the local value of Ji j is, however, not justified in this 
case, because of the existence of long-range correlations 
in the critical region of the XY model. This fact is shown 
in the dependence of the local critical field which does 
not match the spatial profile of the spin-spin interaction 
(Fig. El (a)). In both 1/r 3 and linear trap-XY models, 
B z can be calculated exactly |23, and the result agrees 
with our numerical calculation. 

A better picture of the linear trap-XY model can be 
obtained in terms of Jordan- Wigner fermions |l8| . If one 
neglects long-range terms, which lead to fermion-fermion 
interactions, then a z corresponds to the local density of 
Jordan- Wigner fermions, and the evolution of a z shows 
the emptying of fermionic levels as — B z (which plays the 
role of a chemical potential) is decreased (Fig. 0(b)). 



B. Correlation functions 



Our DMRG calculations show that the phase B z < B z 
is also critical in both 1/r 3 and linear trap-XY models, 
as evidenced in the algebraic decay of Cf J oc l/\i—j\ a . In 
this case, contrary to the Ising model, critical exponents 
are slightly different than in the NN-XY model, in which 
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FIG. 10: Correlation C^j 0+ j between the central ion jo = 26 
and the rest of the chain in the critical region as a function of 
ion separation j. Cfn,jo+j * s we ^ fitted by a power-lay decay 
with the exponents a shown in the figure. 



a = 1/2 (see Fig. 1100 . Note that experiments with linear 
ion traps can detect the algebraic decay of correlation 
functions in XY models, even in the presence of the finite 
size effects induced by variations of Jij . 



V. CONCLUSIONS 

In this work we have discussed in detail a recent pro- 
posal |j| for the realization of quantum spin systems with 
trapped ions under the action of off-resonant standing 
waves. Under certain conditions the coupling between in- 
ternal states and vibrational modes can be written as an 
effective spin-spin interacting Hamiltonian with a resid- 
ual spin-phonon coupling. In this way, the physics of 
quantum criticality can be accessed in experiments with 
ion traps. Our numerical calculations show that: 

(i) In the homogeneous l/r 3 -Ising model, which can 
be realized with ion microtraps there is a quantum 
phase transition with critical field Bf ss 0.83J, and 
same critical properties as the nearest-neighbor 
Ising model. 
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(ii) In linear ion traps, due to the non-constant separa- 
tion of the ions, the spin-spin interaction is inho- 
mogeneous, which leads to the coexistence of differ- 
ent quantum phases. However, critical properties 
can be accessed by measuring each region in the 
trap separately, which only requires partial ion lo- 
cal measurements (10-20 ions). 

(iii) Ion trap Ising models show long-range quan- 
tum correlations that are not present in nearest- 
neighbor models and can be explained by means of 
a spin-wave theory. 

(iv) In the case of the l/r 3 -XY model the critical field 
is shifted to a value sa 0.9. The quantum phase 



diagram can be determined by measuring the effec- 
tive longitudinal magnetization, m z . Besides that, 
experiments with trapped ions can access the alge- 
braic decay of correlation functions in the critical 
phase of XY models. 
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